Interactions and scaling in a disordered two-dimensional metal 
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We show that a non-Fermi hquid state of interacting electrons in two dimensions is stable in the 
presence of disorder and is a perfect conductor, provided the interactions are sufficiently strong. 
Otherwise, the disorder leads to localization as in the case of non-interacting electrons. This con- 
clusion is established by examining the replica field theory in the weak disorder limit, but in the 
presence of arbitrary electron-electron interaction. Thus, a disordered two-dimensional metal is a 
perfect metal, but not a Fermi liquid. 
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A number of recent experiments on the two- 
dimensional electron system in semiconductor devices 
have revealed the existence of a low temperature metal- 
insulator transition as the electron density is varied. In 
the 1980's, similar experiments appeared to confirm the 
prediction that the two-dimensional disordered electron 
system would have no metallic states [||. What distin- 
guishes the new samples from those studied earlier is 
that the electron density is extremely low, so that the 
Coulomb interaction energy is larger than the Fermi en- 
ergy. 

The combined effects of disorder and electron-electron 
interaction were studied by renormalization group (RG) 
methods ||^,^ and there were indications that one effect 
of the interaction might be to stabilize a metallic state in 
two dimensions. However, the nature of the metallic state 
was unknown, and, in any case, the conclusion was not 
very definite since the RG flows went to a strong coupling 
regime beyond the range of validity of the RG equations. 
Therefore, the possibility of a two-dimensional metallic 
state remained an open issue from the theoretical point 
of view. 

A scaling analysis of the recent experiments has been 
made and, based on it, it was argued that any disor- 
dered two-dimensional metal is a perfect metal but not 
likely to be a Fermi liquid. The reason is that if the inter- 
action is turned off, the electrons will localize |^ , and the 
localized state has no resemblance to a Fermi gas. This 
motivated us to consider a two-dimensional non-Fermi 
liquid state and to study the effects of impurities. 

As a minimal specification of a non-Fermi liquid, the 
retarded single particle Green's function must not con- 
tain a quasiparticle pole when analytically continued to 
the lower half plane, but a branch point. This leads to 
a spectral function satisfying the homogeneity relation 
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the non-Fermi liquid. Only the set of exponents yi — 
1, y2 — 1, and y^i = — 1 represents a Fermi liquid for 
which the branch points collapse into simple poles. Here 
the momentum is measured with respect to kp and the 
frequency is measured with respect to the Fermi energy. 

The above spectral function will be assumed to contain 
a kinematic form factor of zero scale dimension, which is 
6{uj^ — Vpk'^). The rationale is as follows: if dissipation 
is due to the decay of an electron coupled to particle- 
hole pairs, then, for a; > 0, a; has to be greater than 
wfI^I- Similarly, because in the ground state of a non- 
Fermi liquid, particles are present both above and below 
the Fermi sea, lu must be less than — for negative 
frequencies. The presence of the 6'-function leads to a 
density of states that vanishes at the Fermi energy. In 
Ref. 0, another choice was made for which the density 
of states remains finite at the Fermi surface. 

A spectral function with non-trivial exponents is nec- 
essary, but it does not fully specify a non-Fermi liquid. 
For example, it does not contain spin-charge separation, 
which requires separate singularities for spin and charge 
excitations. Moreover, new exponents may have to be 
introduced for the scaling of composite operators. It is, 
however, our intention to see what can be learned from 
this minimal specification of a non-Fermi liquid. 

For explicit calculations we use the simple model: 
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in the asymptotically low energy limit, where yi, y2, and 
yA are the exponents defining the universality class of 



where lo^ is a microscopic high frequency scale propor- 
tional to the inverse of the noninteracting density of 
states v. The exponent a has to be positive to satisfy 
the analyticity properties required of the single particle 
Green's function. It is also worth noting that the singu- 
lar part of the spectral function alone will not satisfy any 
sum rules. 

The grassmannian field theory for the localization 
problem was first set up by Efetov, Larkin and 
Khmel'nitskii M and extended by Finkcl'stein p] to the 
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interacting case. We follow Ref. ||]. The random impu- 
rity potential is spatially uncorrelated and has a white- 
noise distribution with zero mean so that {V{r)V{r')) ~ 
(l/27rz^T)5^'''(r — r'), The disorder average is carried out 
on the replicated partition function Zjv, where N is the 
number of replicas. It leads to the action 
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where [tpip) = J2n,ai^ni^n- a^^d Rq 1/4tti^t. The 
Grassmann variables {-^(r), ip"{r)} carry Matsubara in- 
dex n {ujn = n(2n + 1)T); the spin and the replica indices 
are lumped together in a. Here Sq is the action of free 
fermions. It and the electron-electron interaction term 
S'"_o is diagonal in replica indices. Replica mixing takes 
place only in the impurity-induced interaction term. 

If the electron-electron interaction leads to a Fermi liq- 
uid, the impurity induced interaction is relevant by power 
counting. Let us rewrite it as 

RoJ2 [ d'^rdrdr'^" [r,T)^°'{r,T)i^'^ [r,T')ijP [r,T') (4) 

Under the scale transformations k —t Ak and lu Aw, 
the dimension of the fermion variables obtained from the 
free fermion action is [|lO| ~ A~(^/^) so that 

[■(/'"(r, r)] ~ A^-'^/^^ Note that spatial dimensionality, d, 
does not enter into the scale dimensions, because the scal- 
ing is in a direction normal to the Fermi surface; the 
transverse (d — 1) directions act as internal degrees of 
freedom. Then, from Eq. (^, the impurity- induced in- 
teraction scales as R{1) — i?oe', where A"-'^ = e' and 
is thus relevant. Loop corrections show that either the 
system flows to a diffusive metal fixed point (a possibil- 
ity that exists for d > 2) or to the Anderson localized 
phase. The possibility of a pure (ballistic) Fermi liquid 
phase does not exist because disorder is always a relevant 
operator ||ll|. 

In contrast, let us suppose that 5*0 + S^_^ leads to 
a non-Fermi liquid state. Consider the simple model 
in Eq. |. Then it follows that [tfj"{k,uj)] ~ A-(t^), 
and [ij°'{r,T)] - A"^. Consequently, R{1) = i^oe^^"^")'. 
Now the impurity-induced interaction is irrelevant for 
a > 1/2 and relevant for a < 1/2. For a > 1/2, the dis- 
ordered system is described by a non-Fermi liquid fixed 
point |l^]. For a < 1/2, the system has a choice between 
a diffusive metal or a Anderson localized phase. 

To explore this situation further, we rewrite the ac- 
tion in terms of slow variables by isolating the regions of 
small total momentum in the replica-mixed particle-hole 
channel 0. Wc define a matrix B in the replica, spin 

and Matsubara indices by Rq (tpip) = — |Tr£?^(r). The 
Fourier transform of B is 

B:%) ^^=L=yXim^ik + q), (5) 



where f2 is the sample volume. Introducing a hermitian 
matrix field Q{r) to decouple this term we find, using 
a trick due to Schwinger |l^, which is valid even when 
the action is not quadratic in Grassmann variables, that 
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where we have taken the replica limit ^ in a prcf- 
actor, and 
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From Feynman diagrams, one can see that electron- 
electron interaction cannot dress the QB vertices when 
the replica-off-diagonal components of Q are involved. 
However, if Q is diagonal, the vertices can receive cor- 
rections from electron-electron interaction. We shall ne- 
glect the vertex corrections for the replica-diagonal back- 
ground field Q"". This approximation breaks gauge in- 
variance, but has little effect on the localization prob- 
lem, which crucially depends on the off-diagonal Q's. 
With this approximation ||l^, the quantity {TiQB)\ 
can be calculated from the matrix Green's function 
G^'^(r, r'; A), which satisfies the Dyson equation 

G{r,r';X) =G(r-r';0) 

-iX^/2R^ J d'^riG(r-ri;0)g(ri)G(ri,r';A). 

(8) 

The Green's function G(r — r', iuj„i; 0)SmnSaf3 is the fully 
interacting one, but in the absence of the matrix field 
Q{r). It is translationally invariant and diagonal in Mat- 
subara and replica indices. As a result of our approxi- 
mation, G(r — r'; 0) fully determines G(r, r'; A) as a func- 
tional of Q{r). 

Consider the saddle point solution for which Q is a 
constant diagonal matrix, that is, Q{r) = Q{iuJn)6mn5ai3 
M. Then, the self-consistent equation for Q is 



— , i y/2RoG{k,iujn;0) 

- ^ 2^ ^ ^ iQ{lUJ.^)^2IhG{k, iUn] 0) ■ 



(9) 



In the limit {v jr) <C 1, we can write, to leading order. 
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Using the explicit form of the spectral function given ear- 
lier, it is possible to show that 
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where A(a) is a function of order unity. As the next 
iteration, one may substitute this expression for Q in the 
exact saddle point equation to estimate the correction. 
For a > 1/2, it is readily seen that the correction vanishes 
in the limit of zero frequency. For a < 1/2, the correction 
tends to infinity so that Eq. (11) cannot be a correct 
solution. 

For a < 1/2 and not too close to it, the saddle point 
value of Q can be calculated for the model spectral func- 
tion from Eq. (9) without the approximation of Eq. (10). 
For the frequency tending to zero, it is 

1 

Q{iuJn) = B{a) (J^'j sgn(w„) = gsgn(u;„), (12) 

where B(a) is a function of order unity, and we have 
neglected the weak frequency dependence. The saddle 
point solution for a < 1/2 requires careful treatment of 
the cutoff. In the limit a 1/2, the presence of the 9- 
function cutoff in the non-Fermi liquid spectral function 
cannot be ignored, while for a — > it can be. If this 
crossover is treated correctly, one finds that as a — > 1/2, 
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This solution smoothly joins the solution given in Eq. 
(12) at a = 1/2. 

To summarize: In the non-disordered case, the den- 
sity of states for the non-Fermi liquid vanishes at zero 
frequency for all a. In the presence of disorder, this 
quantity is considered as an order parameter, and the 
non- Fermi liquid behavior survives if a > 1/2, that is, 
when the the order parameter vanishes. If, on the other 
hand, a < 1/2, the symmetry is broken, and there is a 
non-zero density of states in the disordered situation. 

Next, we investigate the fluctuation modes about the 
mean-field solution. We consider a small deviation 5Q — 
{Q{r) — Q). It is useful to shift the variable further ||] 
by shifting Q to Q + H, where 

H{iujn) = i\/2-KVT [G"^(fc, iw„; 0) - G-^(k, 0; 0)] /, 

(14) 

where / is the unit matrix in replica indices. Then it is 
not difficult to see that 



S = -T{Q + H)-^ I d'^r TT[6Qf + 2TtH5Q 



-SS, 
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where 5Q = Q — Q — H . The quadratic part of 5S is 
given by 



5S2 = Ro 
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Note that H plays the role of a symmetry breaking field, 
which in the noninteracting case collapses to \/2/rruTijJnI . 
The shift, however, has allowed us to recast the expansion 
in terms of G which is invariant under arbitrary homo- 
geneous unitary rotations of the saddle point. 

The 5Q fluctuation modes will be classified as longi- 
tudinal when mn > and transverse when mn < 0. In 
the Fermi liquid case, the latter are the usual diffusion 
modes. The longitudinal modes are massive as in the 
non-interacting model j^. The zero momentum limit of 
the transverse part of Fm„(0) for a < 1/2 is given by 



G2(fc,0;0) 
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where q is defined in Eq. (^^. One can verify using 
the saddle point equation that -FXi(O) ^ (1/2-Ro)- This 
leads to a complete cancellation between the quadratic 
term — 5 / d''rTr((5(3)^ and 882 in the fluctuation action, 
Eq. (16). The resulting theory for the transverse modes is 
massless, and the remaining quadratic term in the action 
in the long wavelength limit is 



5S2 - 
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where the bare diffusion constant at the scale of the mean 
free path is renormalized by the electron-electron inter- 
action and is (a < 1/2 and not too close) 



D{a) = (l-a) 



TT 
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where D — v^r/d is the usual diffusion constant. In 

this equation we have also normalized the saddle point 
2 

solution such that Q = /. As in the non-interacting 
problem, the resulting theory can now be cast in the lan- 
guage of a matrix nonlinear cr-model. It is well-known 
that all states are localized for d < 2. The only differ- 
ence is that the bare diffusion constant D{a) is enhanced 
over the non- interacting value D. 

As in the corresponding non-interacting problem, the 
mass term cancels in the transverse channel and the dif- 
fusion propagator diverges as fc~^ at zero frequency, for 
a < 1/2. These are properties of the Goldstone modes 
in the transverse replica channel and follow as long as 
the continuous symmetry in the Q-field theory is broken, 
which is indeed the case for a < 1/2; they do not depend 
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critically on the assumptions made in the longitudinal 
channel. Thus, the neglect of the vertex correction in 
the longitudinal channel does not spoil this picture. 

For a > 1/2, Q{iujn) tends to zero as the frequency 
tends to zero, and the transverse part of the mass term 
does not cancel against the transverse part of SS2 ■ Our 
analysis shows that impurity interactions are irrelevant 
for a > 1/2, and the theory corresponds to the non- Fermi 
liquid fixed point Hamiltonian, minimally specified by 
the scaling property of the one-particle spectral function. 
It is preferable then to return to the original fermion 
language. Such a system should be a perfect metal. 

It is important to make a number of qualifying re- 
marks. 

(1) For a > 1/2, any attempt to expand in powers of 
the exponent a will lead to numerous logarithmic diver- 
gences, and the theory cannot be cast in the language 
of a Fermi liquid theory. For a < 1/2, the electrons are 
localized in d < 2, with the interacting bare diffusion 
constant enhanced over the non-interacting value. Thus, 
once again, the language of Fermi liquid theory is not 
meaningful. 

(2) The model spectral function considered here does 
not contain spin charge separation. From a simple model 
reflecting spin charge separation [Q, it is possible to see 
that our qualitative conclusions regarding broken sym- 
metry remain unchanged. It is not apparent, however, 
that that is the correct model for spin-charge separation. 
In fact, the entropy producing process in which an ac- 
celerated electron breaks up into a spin and a charge 
excitation is missing iQ. It is also possible that spin- 
charge separation will affect the behavior in a magnetic 
field. For example, in a two-dimensional system, a par- 
allel magnetic field can couple to the spinon excitations, 
producing a gap in the spectrum, thereby destroying the 
non-Fermi liquid state. This, in turn, will lead to a lo- 
calized insulating state as seen in experiments |l5|l . 

(3) The results obtained here are only superficially sim- 
ilar to those of the one-dimensional problem |l^ . In one- 
dimensional theories, the anomalous dimension of the 
density-density correlation function cannot be simply ob- 
tained from the anomalous dimension of the one-particle 
Green's function. As remarked earlier, the operator prod- 
uct expansion of composite operators contain further sin- 
gularities. Such possibilities were purposely avoided to 
reduce the complexity of the theory. Nonetheless, we 
have demonstrated that the breaking of the single most 
important assumption of the Fermi liquid theory, namely 
the existence of the quasiparticles, leads to dramatically 
different behavior with respect to the localization prop- 
erties of the interacting electron system. Moreover, it 
appears that there are experimentally relevant regimes 
in two dimensions, which show no hint of instabilities 
that could arise from 2fci? vertex corrections. Thus, the 
neglect of vertex corrections in two dimensions may have 
a well-defined range of validity. 
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